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Abstract. In this paper, we compute the Laplace transform of occupation times (of the nega- 
tive half-Une) of spectrally negative Levy processes. Our results are extensions of known results 
for standard Brownian motion and jump-diffusion processes. The results are expressed in terms 
of the so-called scale functions of the spectrally negative Levy process and its Laplace exponent. 
Applications to insurance risk models are also presented. 



1. Introduction and results 

Let X = (Xf )t>o be a spectrally negative Levy process, that is a Levy process with no positive 
jumps. The law of X such that Xq = x is denoted by P^; and the corresponding expectation by 
E3.. We write P and E when x = 0. As the Levy process X has no positive jumps, its Laplace 
transform exists and is given by 

E 



for 9,t > 0, where 

m = iO + \y''Q' + y"° (e'^ - 1 - 0^i(_i,o)(^)) n(dz), 

for 7 G M and cr > 0. Also, 11 is a u-finite measure on (—00, 0) such that 

/■° 

/ (1 A z2)n(dz) < cx). 

J —00 

The measure 11 is called the Levy measure of X, while (7, o", 11) is referred to as the Levy triplet 
of X. Note that E \X{\ = '(A'(0+). 

For an arbitrary spectrally negative Levy process, the Laplace exponent ^ is strictly convex 
and lime_).oo ^(0) = 00. Thus, there exists a function ^ [0, 00) — )• [0,oo) defined by ^"(0) = 
sup{^ > I ^(^) = 0} (its right-inverse) and such that 

V;($(0))=0, 0>O. 

We first examine the total occupation time of the negative half-line (— oo,0). 

Theorem 1. If'tp'{0+) > 0, then, for A > 0, 



(1) E re"^/o°°l{x.<o}dsl ^ tP'{0+) 



,$_(A) 
A 
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where <I>(A)/A is to be understood in the limiting sense when A = 0. 

We now recall the definition of the g-scale function W^'^\ For q > 0, the (^-scale function of 
the process X is defined as the function with Laplace transform on [0, oo) given by 

/■oo 1 

/ e-^'W^''\z)dz = — , for 9 > <^{q), 

Jo w{G) - q 

and such that W^'^\x) = for x < 0. This function is unique, continuous, positive and strictly 
increasing. We write W = W^^^ when q = 0. We have that W^'^^ is differentiable except for at 
most countably many points; see pTBj . Moreover, H^('J) is continuously differentiable if X has 
paths of unbounded variation or if the tail of the Levy measure, i.e., the function x i— n(— oo, —x) 
on (—00,0), is continuous. Further, W^'^^ is twice continuously differentiable on (0, oo) if o" > 0. 
For more details on the smoothness properties of the g-scale function, see [4J. We will also use 

the functions | W^'^^ g > o} and {Z^'^^;q> O} defined by 

W^''\x)= rW^''\z)dz. 
Jo 

and 

Z(5)(x) = l + qW''''\x). 
We can now state the following corollary to Theorem [TJ 

Corollary 1. //t/;'(0+) > 0, then, /or A > and x >0, 

/•oo 

= ^'(0+)$(A) / e-'^^^^'W{x + z)dz, 
Jo 

Note that when x = 0, by the definition of scale function, we have 

Jo H'^W A 

therefore recovering Theorem [T] as a special case. 

These two results generalize the work done in |21| where the sum of a compound Poisson 
process and a Brownian motion is analyzed (see, e.g.. Equation (4.9) in that paper). Using ruin 
theory terminology, they study the duration of negative surplus, also called the time in red, 
in such an insurance risk model; their work is itself an extension of ^ in the pure compound 
Poisson case. 

We now examine the occupation time of (— oo,0) until a negative level —b is crossed for the 
first time. Let be the first passage time below —b of X: 

tZ^ = mi{t >0: Xt<-b}. 

Theorem 2. If'tp'{0+) > 0, then, for A > 0, 

nO+) + 4^ + /°- 4'\x) Q U{dx - y)dy 



(2) 



E 



{Xs<0} 



ds 
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where 




Af\h) = zW(6)TyW(6) - A(TyW(6))2, 



and 




As a special case, we recover the corresponding result for standard Brownian motion; see, e.g.. 



Our proofs use fluctuation identities for spectrally negative Levy processes and, as a con- 
sequence, our results are expressed in terms of the so-called scale functions of the spectrally 
negative Levy process (see, e.g., [T| [2]) and its Laplace exponent. Only elementary arguments 
are needed. To the authors' knowledge, the literature seems rather scarce on the relationship be- 
tween scale functions and certain occupation times of a general spectrally negative Levy process. 
However, similar results can be found in |12j . 

The current work has been partly motivated by the study of an insurance risk model with 
implementation delays. Insurance risk models use stochastic processes to describe the surplus 
of an insurance company. In risk models of a Parisian nature, an implementation delay in the 
recognition of an insurer's capital insufficiency is applied. More precisely, it is assumed that ruin 
occurs as soon as an excursion below a critical level is longer than a deterministic time; such 
models have been studied only very recently in |6l [71 [T7] and the idea has been borrowed from 
finance and Parisian barrier options (see [5]). Of more interest in our context is the work by the 
same authors: in |17| . instead of a deterministic delay, an exponentially distributed grace period 
is used in the definition of the Parisian ruin. It turns out that the probability of ruin in this 
model is strongly related to the occupation time of the underlying process. The reader is invited 
to consult Section [6.21 to obtain further details on this connection. 

The rest of the paper is organized as follows. In the next Section, we recall the relevant 
notions and results on scale functions and fluctuation identities. Then, in Section 3 and 4, we 
prove Theorem [2] and Theorem [1] respectively. Section 5 presents applications to the case of a 
Brownian motion with drift and insurance risk models. 



mm- 



2. Scale functions and fluctuation identities 



We recall some of the properties of the g-scale function T^^'?) and its use in fluctuation theory. 
Let d = 7 - zn(dz). The initial values of and VF^?)' are known to be 




and 



'2/<T^ when cr > 0, 

VF('?)'(0+) = I (n(-oo, 0) + q) /(f when = and n(-oo, 0) < oo 



oo otherwise. 



Now, define 



~ 



inf{t > 0: < 0}, 
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and, for a > 0, 

r+ = mf{t >0: Xt>a}, 
with the convention inf = oo. It is well known (see, e.g., |15| ) that, for x < a, 



-QTa 



E. 



Wii){a) 



V(5)(a)' 



and 



-4>(g)(a-iO 



If 'ip'iO+) > 0, then P^; {t+ < oo} = 1 and therefore 
Also, we have that 

2 

(3) {v < oo,X^- G dz} = y [H^'(:r) - ^{0)W{x)] 6o{dz) 

+ ^ U{dz - y) ^e-'^^^'^yW{x) - W{x - y)} dy, 

where 5q is the Dirac measure at 0. The first term of this measure corresponds to the case when 
X^- = 0, a behaviour called creeping. Note that we write to mean J^^^^y similarly, we will 

to mean L ni , etc. 



{tq < oo} = 1 
n(dz - y) {W{x) - W{x - y)} dy. 



write to mean /(_^ q), 

Using the distribution in Equation (|3]), together with the fact that 
ip'{0+)W{x) if ^'(0+) > (in which case ^>(0) = 0), we obtain 

2 

(4) 1 = ijj'{0+)W{x) + —W'{x) + 

2 J-oo J0+ 

For more details on spectrally negative Levy processes and fiuctuation identities, the reader 
is referred to [15j . Further information, examples and numerical techniques related to the com- 
putation of scale functions can be found in [ll|9l [101 ESI [20]. 

3. Proof of Theorem [2] 

The main idea of the proof consists in defining a quantity underestimating and overestimating 
the occupation time 

I{x,<o}ds. 



This respectively leads to an upper and a lower bound to its Laplace transform. Subsequently, 
by taking an appropriate limit, we show that the two bounds converge to the expression on the 
right-hand side of ([2]). 

First, we provide a lower bound to this Laplace transform by overestimating the occupation 
time. To this end, we consider a clock which starts at time and stops when level e is attained 
or when level —b is crossed. Then, if level e was attained first, every time we go below 0, we 
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re-start the clock and subsequently stop it when we get back to e (without going below —6); let 
L^''' be the Laplace transform of this overestimating quantity of the occupation time when the 
process X sits at level e at time 0. Hence, by the strong Markov property of X, we have 



E 



> E 



Using the strong Markov property and the spatial homogeneity of X, we get 



L 



,{tq =oo} 



+ y ^Pe {tq- < 00,X^- e dx} {e. 



E., 



e'^"' ;r+ <T_, 



Note that when x < —b, 



and 



Consequently, 



E, 



E, 



L 



F,{t- =oo} + /°+ P, {t- < oo, X^- e dx} E,. 




_ Lf + 


l-/°iPe{ro-<oo,X-Gdx}E, 




^3 



where, using some of the fluctuation identities in Section [2j 



Ll'" = ^lJ'{0+)W{e) + ^W'{e 



L 



and. 



zW(x + 6) - z(^)(6 + e) 



W^^\x + b) 



2 ^'VW(6 + e) 



-oo WW{b + e) Jo 



WW{b + e) 

/•oo 

/ U{dx-y){Wie)-W{e-y)}dy, 



Uidx-y){Wie)-Wie-y)}dy. 



Using Equation (jH, we get 



L 



W'je) 
YwW{b + e) 



+ 



r 

J —ex. 



1 



W'^^\x + b) 



W^ib + e) 



U{dx-y){Wie)-W{e-y)}dy. 



Alternatively, we now develop a scheme to underestimate the occupation time in question. 
Every time we go below — e (— e > —6), we start the clock and stop it when we get back to 
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(without going below —b); let t/*^'^ be the Laplace transform of this underestimating quantity of 
the occupation time when Xq = 0. Hence, by the strong Markov property of X, we also have 



E 



< u 



As above, using the strong Markov property and the spatial homogeneity of X, we can write 



+ 



e{ro = 00} 



< oo,X - G dx N E^. 



}{■ 



—Xt , , — ^ -\- 



+ E^ 



and then 



where 



U' 



M _ ^1 + ^2 



Z^^\h - e) - Z^^\b) 



U. 



and 



0- 



0+ 



As for Lg^, using Equation (jl]), we get 



Il{dx-y){W{e)-W{e-y)]dy, 
Il{dx-y){W{e)-W{t-y)]dy. 





+ / 1 



W^^\b)-W^^\b-e) 

W^^\x + b-e 
WW(h) 



Il{Ax-y){W{e)-W{e-y)]dy. 



3.1. Proof if cr > 0. First, we assume that X has a Brownian component, that is cr > 0. In 
this case, VF(0) = and W^(^)'(0+) < 00. Then, 



/(0+)H^'(0+) + 



W'{<d+) 

2 



lim — ^ = lim 

e-i-O e e-^O 



P^'(0+) 



zW(x + 6)-ZW(6 + e) 



n(dx — y) 
0+ I f 



n(dx - y)dy 



0+ 



OCCUPATION TIMES 



and, similarly, 

lim = V (0+) lim — — + — lim ' , 



+ lim 



n(dx - y) 



0+ 



W{e) W{e - y) 



dy 



^'(0+)H/'(0+) + -W"(0+)-^ 



+ W'{0+) / 

t/ — ( 



WW{b) 



U{dx - y)dy. 



0+ 



The result follows. 



3.2. Proof if cr = 0. If X does not have a Brownian component, we adapt the previous method 
as follows: 

re,b 

V'(0+), for all e > 0, 



lim ^ = lim 



lim 

e->0 



zW(x + 6)-ZW(6 + e 



zW(x + 6)-.zW(6 + e 



VW(6 + e) 



zW(x + 6)-ZW(6) 



oo 

0+ 
oo 

0+ 



n(dx - y) 

Il{dx - y)dy 



W{e) W{e - y) 
W{e) W{e) 



dy 



and, similarly, 

lim — ^ = V'(0+) + hm 

= V'(o+) + 

The result follows. 



WW{b) 

^ W'^^^x + b) 
~ WW{b + e) 

W^^^x + b) ] r 
WWib) Jo+ 



poo 

/ Il{dx - y)dy 



0+ 



U{dx - y) 



W{e) W{e-y) 
W{e) W{e) 



dy 



Ti{dx - y)dy. 



3.3. Conclusion of the proof. In all cases, i.e., when X has or does not have a Brownian 
component, the limiting results for V^'^ can be obtained in a similar fashion. In fact, it can 
easily be proved that, for each i = 1, 2, 3, 



hm — 

L^. 



e,b 



and then 



e->o wW{b + e) 
The details are left to the reader. 



lim U''K 
e->0 



8 landriault, renaud and zhou 

4. Proof of Theorem [T] 

We use the same general idea as in the proof of Theorem [2] Define L*^ = L^'°°. In this case, 
using the strong Markov property twice, we obtain 

rO+ 



= F,{tq = oo} + L'E e'^''^ J ^^t^ < oo, X^- G dxj 
= V''(0+)VF(e) + L^e-*(^)^ / I < oo, X^- G dx \ e*^^)^' 

J — oo J 



Since, for r > 0, 
(5) 

we can write 

fO 



r rX _ _ 
e ^0 ; Tq < oo 



Jo 



W{x), 



j {tq" < OO, X^- e dx} e*(^)^ = E, 



<I>(A)X _ 

e ^0 ; Tq < oo 



It is easily shown that 



lim , ^ ' = 0. 



Indeed, if VF(0+) > this is trivial and if M^(0+) = one can use L'Hopital's rule. Consequently, 



lim L*^ = lim 



^'(0+)M^(e) 



1 



1 - A /o e-^W^W{x)dx - ^e-*W^ W(e) 



Similarly, if we define = U^'^, one can show that 



lim U' = ij'{0+)^^ 

e->-0 A 



and the result immediately follows. 



5. Proof of Corollary [T] 

Using the strong Markov property of X twice, using some of the fluctuation identities in 
Section [2] and then using Theorem [H we have 



-^/o°°I{x^<o}ds 



^^{tq- =oo}+E^. 



i;'iO+)Wix) + ij'{0+) 



— oo 




Tf) < oo,X - G dz E 



} 



e -^^0 ; ro+ < oo 



i;'{0+)W{x) + ^^'(0+)^ I Px- [tq < oo,X^- G dz} e*^^ 



V''(0+)VF(x) + V''(0+)^E, 



e ^0 ; Tq < oo 
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Using the result in Equation ([5]) once more, we get 



Ex 



-■^/o°°1{Jfs<o}<is 



)} 



V,'(0+)i^|e*W- (l-A 
= V'(0+)^ e-*W"t^(x + z)dz 

where in the last step a change of variables and an integration by parts were undertaken. The 
result follows. 



6. Applications 

6.1. Brownian motion with drift. For Brownian motion with drift, i.e., if X is of the form 
Xt = mt + aBt with m = ^^'(0+) > and cr > 0, where i? is a standard Brownian motion, then 
n = and, using Theorem [2l we have 



E 



mWW{b) + {a^/2)WW'{b) 
where, as one can easily verify from the definition. 



^/m? + 2A(t2 
for X > 0. 

Then, for standard Brownian motion, in which case m = '(/''(0+) = and o" = 1, we have, 
from Theorem [21 that 



E 



where 



and then 



— sinh(xv 2A) 



= cosh(xV2A). 



Recalling that cosh^(x) — sinh^(x) = 1, we get that 



E 



2cosh^(6V2A) - A(2/A) sinh^(6V2A) 



2cosh(6V2A) 



cosh(6V2A) 

therefore recovering Proposition 4.12 of [l3] (see also |llj). 
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6.2. Insurance risk models. Classical insurance risk models describe the surplus process of 
an insurance company using a compound Poisson process or a Brownian motion with drift, i.e., 
special cases of spectrally negative Levy processes. In those models, it is usually assumed that 
the net profit condition holds; this condition ensures that ruin will not occur almost surely. For a 
Levy insurance risk process, i.e., a spectrally negative Levy processw with non-monotone paths, 
it amounts to E[Xi] = iIj'{0+) > and it is also equivalent to liuit^oo Xt = oo almost surely. 
An interpretation of Levy insurance risk models for the surplus modelling of large insurance 
companies is for instance given in ^14j . 

Classically, the probability of ruin has been the most studied risk measure to evaluate the 
quality of a company. More recently, the analysis of the duration of the negative surplus, in other 
words the occupation time of the negative half-line, has also been considered in the compound 
Poisson case ^ and then in a jump-diffusion model [21j; see also |3]. Recall that Theorem [1] 
extends Equation (4.9) in |21| where the sum of a compound Poisson process and a Brownian 
motion is considered. 

We now want to provide another link between Theorem [1] and insurance risk models. In ^17j . a 
new definition of the time to ruin is proposed. In that paper, each excursion of the surplus process 
X below is accompanied by an independent copy of an independent (of X) and exponentially 
distributed random variable e^ with mean l/d; we will refer to it as the implementation clock. 
If the duration of a given excursion below is less than its associated implementation clock, 
then ruin does not occur. More precisely, we assume that ruin occurs at the first time Td that 
an implementation clock rings before the end of its corresponding excursion below 0. It is worth 
pointing out that the time to ruin is easily defined when the Levy insurance risk processes X 
has sample paths of bounded variation. 

Therefore, in |17| . the case of a surplus process of bounded variation is considered. In that 
model, one can show that the probability of ruin in this Parisian risk model with exponential 
implementation delays can be expressed as follows: 

¥{Td < oo} = l-E\e'~'^^o"°^l^s<o}'^A ^ 

if the exponential clock has mean 1/d. Indeed, for each excursion of length T, the probability 
to survive, i.e., the probability that the exponential random variable associated with it is larger 
than T, is equal to exp{—dT}. Using the independence assumption between the clocks and 
summing up over all the excursions (there are countably many of them), we get 

Using Theorem [H we recover the corresponding expression in |17) , that is 

P{r,<oo} = l-V'(0+)^, 

when the net profit condition tp'(0+) > is verified. 

More generally, using Ito's excursion theory for spectrally negative Levy processes, this Parisian 
risk model with exponential implementation delays can also be defined when the underlying sur- 
plus process has paths of unbounded variation. It suffices to mark the Poisson point process of 
excursions away from zero (see [19] for a definition of the corresponding excursion process) with 
independent copies of the generic random variable e^, similar to the proof of Theorem 6.16 in 
[15) : for an excursion away from zero starting above zero, this time spent above zero is simply 



OCCUPATION TIMES 



11 



ignored. As a consequence, Theorem [T] provides a generalization of the probability of ruin in a 
general Levy insurance risk model with exponential implementation delays. 
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